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We discuss the existence of monodromies associated with the singular points of the eigenvalue 
problem for the Rabi model. The complete control of the full monodromy data requires the taming of 
the Stokes phenomenon associated with the unique irregular singular point. The monodromy data, 
in particular the composite monodromy, are written in terms of the parameters of the model via 
the isomonodromy method and the tau-function of the Painleve V. These data provide a systematic 
way to obtain the quantized spectrum of the Rabi model. 
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INTRODUCTION 


The Rabi model 01 describes the interaction of a sin¬ 
gle electromagnetic mode - a simple harmonic oscillator 
- with matter - a two-level quantum system. It is a quite 
simple model yet it has an interesting and rich spectrum. 
Recently, it has attracted vividly attention due to ex¬ 
perimental and mathematical reasons. From applied and 
experimental physics side, there have emerged interesting 
applications in quantum optics and quantum computa¬ 
tion because of good prospects of experimental realiza¬ 
tion using Josephson junctions, traped ions and others 
(see references in i)- From the mathematical side, for 
a long time it has been a challenge to prove its exactly 
solvability. And so, finally, as recently as 2011, Braak [3| 
solved the model in a Bargmann (coherent state) repre¬ 
sentation obtaining in a systematic way its regular spec¬ 
trum as zeroes of a transcendental function. The excep¬ 
tional part of the spectrum were already known since the 
late 1970’s 0, but it can also be obtained in the frame¬ 
work proposed by Braak ! 5j. 

This fact has opened up a whole new set of interest¬ 
ing problems in mathematical physics. One issue regards 
the notion of integrability of the Rabi model. At face 
value, the Rabi model could be a first instance of an 
exactly solvable yet not integrable model in mathemat¬ 
ical physics. A controversy ensues since the proper def¬ 
inition of integrability in quantum physics - or even in 
mathematics - seems not to be clear cut (ft [7]. Braak 
claimed that the Rabi model is integrable in a new quan¬ 
tum integrability criterion coined by himself. A system 
is dubbed Braak integrable if there are f = f c + fd quan¬ 
tum numbers classifying the eigenstates uniquely, where 
fc, fd stand for the number of continuous and discrete 
degrees of freedom (d.o.f.), respectively. For the Rabi 
model, this number is f c = fd = 1 , one harmonic oscilla¬ 
tor and one two-level system, so that / = 2. These two 
quantum numbers arises from the Z 2 parity symmetry of 
the Hamiltonian. Even in the lack of a better mathemat¬ 


ical formulation of this criterion - still a pending project 
- it served as the guiding principle to the solution of this 
outstanding problem. 

Controversies aside, it would be desirable to frame 
this discussion in a more conservative setup. Therefore, 
Batchelor and Zhou [j] raised the issue of whether the 
Rabi model is Yang-Baxter integrable (YBI). They suc¬ 
ceed to show YBI for two special points in the space of 
parameters of the Rabi model. For a generic point in this 
moduli, YBI is still an open question. 

Another important issue arising from the Braak’s work 
is that of finding an universal method that applies to 
a wide variety of models involving coupling of a boson 
mode with a two-level system in the Bargmann repre¬ 
sentation. Such a program has been developed by Ma- 
ciejewski et al. f§, 9j • They have used a method framed 
in terms of the Wronskians, that is, a 2 x 2 matrix con¬ 
taining both the wave-function and its derivative in the 
neighbourhood of a singular point. 


The present work advances some enlightenment in the 
direction of proving YBI for a generic point in the moduli 
of parameters of the Rabi model. The departing point is 
the Bargmann representation of the Rabi model. In this 
representation, it can be easily shown 0 that the Rabi 
model is described by a confluent Heun equation. Here 
we use the known fact 01 [ll| that the monodromy data 
of these equations can be cast in terms of Painleve V 
transcendent tau-function 12] via isomonodromy equa¬ 
tions. The global properties, relevant for the problem of 
the spectrum and for the YBI, are encoded in the notion 
of composite monodromies. We then present the com¬ 
posite monodromy parameter of the Rabi model. Finally 
we discuss in general terms how one could use this com¬ 
posite monodromy parameter to obtain the Rabi model’s 
spectrum. 


The novelty of our work consists in the presentation of 
the monodromies associated with the singular points of 
the ODE arising from the eigenvalue problem of the Rabi 
model. We then discuss the relevance of the Stokes phe- 
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nomenon in order to have a complete monodromy data 
set. We conjecture that it is the emergence of the Stokes 
phenomenon and the need of extra parameters in the 
monodromy data set that rendered extra difficulties in 
the full demonstration of the Yang-Baxter integrability 
of the Rabi model. 

This work is organized as follows: we first write the 
Rabi model as a standard Fuchsian system, and then 
discuss the mododromies around the singular points. A 
special situation happens for the monodromy around the 
unique irregular singular point, the point at infinity, giv¬ 
ing rise to the Stokes phenomenon, which we discuss in 
detail. As an outcome we obtain the general group re¬ 
lation for the monodromies and how it is related with 
Yang-Baxter equations. We next discuss the isornon- 
odromy method aiming at writing the composite mon¬ 
odromy parameter in terms of the monodromy parame¬ 
ter at the irregular point and the stokes parameters. On 
the other hand, the existence of monodromy matrices for 
our original system are obtained from the tau-function 
of the Painleve V. We finally obtain the composite mon¬ 
odromy parameter in terms of the parameters of the Rabi 
model which provides in a systematic way the quantized 
spectrum of the Rabi model. 


RABI AND ITS MONODROMIES 


The quantum Rabi model (QRM) is described by the 
Hamiltonian 

Hr = a^a + A a z + ga x (a 1 ^ + a) , (1) 

where the boson mode is described by [a, ad] = 1, the 
fermion mode by the Pauli matrices, A is the level sep¬ 
aration of the fermion mode and g is the boson-fermion 
coupling. 

The QRM can be written in terms of two copies of 
Jaynes-Cummings model (JCM), each with its appropri¬ 
ate chirality. Indeed, the chiral JCM Hamiltonian reads 

Hjc = o^a + A u z + g (cr + a + a~a^ , (2) 

and the anti-chiral one reads 

H-jq = a 1 a + A <j z + g (er _ a + cr + a^) , (3) 

so that 

Hr = ^ ( Hj C + Hjc) = + Act 2 + ga x (a [ + a), (4) 

where we have used a x = ( a + + a~)/2. It is important 
to notice that [Hjc, Hjc\ ^ 0- 
Consider the Ansatz !$, Q 

Wa t ))=/ 1 (at ) |0)|+ ) +/ 2( at ) |0)|->, (5) 


where the harmonic oscillator ground state is defined by 
a|0) = 0, cH|±) = ±|±), and /*, i = 1,2, are analytic 
functions of a'. We can now use Bargmann’s prescription 

a 1 i —> w, a H > d w , (6) 


so that [a, a']f(w) = f(w). Substituting the Ansatz into 
the stationary Sclrroedinger equation, or the eigenvalue 
equation, Hr\i/j) = E\ip), we obtain, after setting f± = 

fi ± h, 


E-gw A 

o w f+ = -;-/+-—/- 


dwf- = - 


w + g 
A 
w- g~ 


f + 0LtJE f 

w - g 


( 7 ) 


Eliminating say f + in terms of /_, this results in a second 
order linear differential equation for f + (and /_) which 
could be brought to a confluent Heun equation [5(. Let 
us define 


2 = -2 g(w + g), 






( 8 ) 


( 12 ) 

where /^_ ’ ; are the two linearly independent solutions 
of the system above. The fundamental matrix c f > ( 2 ) is 
then invertible and unique up to right multiplication of 
a constant matrix. With this change of variables, we can 
bring the model to a standard Fuchsian form: 


-$ 1 = —(T-? 

dz 2 3 


1 . 1 . 

— Ao H-- A*, 

z z — t 


with t = —4 g 2 and 


^3 



e1 9 2 )' 



( 9 ) 


( 10 ) 


The system © has two regular singular points at 
Zi = 0,t, i = 0, t and an irregular singular point at 
Zoo = oo with Poincare index 1. The analytical structure 
of the system near the regular singular point is charac¬ 
terized by the monodromy matrices and M t , defined 
as the effect of an analytical continuation around the cor¬ 
responding singular point, 

$((z - Zi )e 2ni + Zi ) = $(z)Mi, i = 0,t. (11) 

One notes that, since any two set of solutions 4>(z) are 
related by right multiplication, the monodromy matrices 
are defined up to an overall conjugation. Moreover, one 
can choose initial conditions of © so that, near a regular 
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singular point Zj, one has 


$(z)| 


Z~Zi 


({z-Zi) a t 0 \ 

V 0 (z - Zi) a * ) ' 


( 12 ) 


Therefore, one can see that generically the nronodromy 
matrix can be written as 


Mi = gi 




(13) 


where <ji G SL(2, C) are called the connection matrices. 
They are also defined up to left multiplication. One also 
notes that, for algebraic purposes, only the difference 
Oi = — ce~) is important. Overall shift of the coeffi¬ 

cients af can be obtained by an s-transformation of the 
solutions: f±(z) -1 (z - Zi) a f±(z). 

The system at the irregular singular point at z = oo 
is slightly more complicated, due to the Stokes phe¬ 
nomenon. In order to describe it, let us start by not¬ 
ing that close to z = oo the solutions are of the form 
f±’ 2 ^ ~ e ±£ / 2 , and the Frobenius series obtained at this 
point is only formal: its convergence radius is zero. One 
can see that generically near infinity the system has the 
form 


where the numbers Sk are called Stokes parameters. By 
the identification of Sj +2 with 2 - 7 T rotations, one then 
defines the nronodromy at infinity at sector Sj by 

Moo\ Sj = SjS j+1 e™ 9 °°° 3 , (19) 

which is the one satisfying the free group relation with 
the other two monodromy matrices: 


MooMt.Mo = 1. (20) 

Note that once we settle in a sector Sj, say j = 1, then 
knowledge of only two consecutive Stokes parameters si 
and S 2 are sufficient to reconstruct the whole series of 
Stokes matrices Sj. 

The outcome of the above analysis is that the whole 
set of parameters 6 = {do, $oo, si, S 2 } is sufficient to 
determine the monodromy matrices up to an overall con¬ 
jugation. This set is thus called the monodromy data. 

The existence of the monodromy matrices provides an 
explicit representation of the 3-braid group in fact, the 
permutation group S 3 - acting on M,; as 

aij(Mi) = MjMiMj\ 

( 21 ) 

a ij {Mj ) = (MjMJMjiMjMj)- 1 . 


d<S>^ 1 

dz 2 3 


An 


+ 0(z~ 2 ). (14) 


From the coefficient A aQ = — (Ao + A t ) of the z _1 term 
one can define the naive monodromy at z = 00 , given by 
the difference of the eigenvalues 9oo = The 

epithet “naive” comes in because the monodromy struc¬ 
ture around z = 00 also depends on the first constant 
term. To describe this structure, we follow 1(1 12. 131 
and define the sectors of the complex plane: 


Sj = {z G C | (2 j - 5)| < argz < (2 j - 1)|}, (15) 


j = 1, 2,.... On each Sj we have the following asymptotic 
behavior for the solutions of the system ©: 


$(z)|. =Gj(z- 1 )exp(±za 3 )z-5 e °°™, (16) 


where G(j(z -1 ) = 1 + 0 (z -1 ) is analytic near z = 00 . 
The Stokes phenomenon relates the solutions satisfying 
(fTUD between different sectors Sj , 


$ 7 + 1 0) = ®j(z)Sj, (17) 

where Sk are the Stokes matrices. Now, $j(e 2m z) = 
&j+ 2 (z)e- m6 °° a3 - they are defined in the same domain 
-, so we have that S j+2 = e m9,x>173 Sje m9a ° a3 . Therefore 
one can choose a basis where 


These generators satisfy 


&ij O (Jjk O — &ik 0 &kj 0 ®ji") 


( 22 ) 


which is known as the Yang-Baxter relation [14j. The 
existence of the monodromy matrices then assures that 
the Rabi model is integrable in the algebraic sense. 


THE ISOMONODROMY METHOD 


The Riemann-Hilbert problem consists in finding a 
Fuchsian system with a prescribed set of monodromies. 
Our problem is the inverse one. In order to solve such in¬ 
verse Riemann-Hilbert problem, we will first notice that 
there are many different families of Aj’s which give the 
same monodromy. Some of them are trivially related by 
overall conjugation. But even so there is still a family of 
non-trivial set of Fuchsian systems parametrized by the 
position of an extra singular point t. 


This family was first described by Schlesinger - see [11 
3 EE El for reviews - but it is more easily understood 
in terms of flat holomorphic connections 17]. Suppose 
we set 


A{z,t) = -03 + 


Ap(t) 

z 


At(t) 
z — t 


d<b{z,t) 

dT~ 


(23) 


S 2J = 


1 «2 j 
0 1 


S'2j +1 ^ 


1 0 

S2j+1 1 


(18) 


as the “z-component” of a flat connection. It is straight¬ 
forward to see that, if we consider the “i-component” 
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relation between the parameters: 
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B(z,t) = 


At{t) 

z — t' 


(24) 


- 1 Ot - 1 


A — £ 


n- 


C 0 


Ct 

A — £ 


(30) 


then F = dtA — d z B + [ A , B) = 0 if the Aj(£) satisfy 


dAt 

dt 


dA ° _ 1 M A 1 

— - -[A t ,A 0 ], 

= --^[At,A 0 \ - -[A t , 0 - 3 ]. 


(25) 


This system are called the Schlesinger equations. Since 
the “field strenght” F vanishes, then the monodromy 
data of the Fuchsian system |9]) will be independent of 
t if A 0 (£) and A t (£) satisfy the equations (|25|) . The ma¬ 
trices Aq and A t can be thought of as a Lax pair for the 
isomonodromy flow. 


Despite being seemingly more complicated, the 
Schlesinger equations (E5l) have a Hamiltonian structure. 
The most direct way to illustrate it is to consider the 
EDO associated with the generic Fuchsian system ([23]) . 
Let us choose a gauge such that 

TrA oo =0 oo . (26) 


This relation means that the change of t has to be ac¬ 
companied by a change in A and /j, so that the relation is 
maintained. 


The Schlesinger system (USD, when written in these pa¬ 
rameters, yield the Painleve V equation for the following 
function of the entries of the Af. 

_ (Ao)n(At)i2 _ + Ot — Oqo — (2 p — 1)(A — t) 

^ (At)n(Ao)i2 6*o + dt — #00 — (2/i — 1)A 


The Painleve V is part of the family of second order 
differential equations with rational coefficients and the 
Painleve property: all the branch points of the solutions 
are fixed, determined by the ODE itself 11 ,0. These 
equations define new special functions, and the Painleve 
V in particular has been useful to compute correlation 
functions of strongly coupled Bosonic systems 0 , dis¬ 
tribution functions of random matrix theory, certain lim¬ 
its of conformal blocks and the XY model - see [20] for 
a (not exaustive) list of applications. It has also been 
shown to give exact analytic expressions for the scatter¬ 
ing of massless fields in black hole backgrounds 21, 22| . 


Consider the off-diagonal term Ai 2 of (E51) . It is of the 
form 


A i 2 ( 2 ) 


k{z — A) 
2(2 — t) ' 


(27) 


where k, A are linear functions of (Ao)i 2 and (Af)i 2 . Now, 
by writing the solution as 


We follow [l3| and define the tau-function: 

logT(£, 0) = — i Tr cr 3 A t — i Tr A 0 A t (31) 

which satisfies a third order non-linear ODE - the so- 
called cr-form of the Painleve equations and it is defined 
up to a multiplicative constant. 


$(2) 


(f? w /f w) 

U a) « 


(28) 


one can check that the elements of the first row f+’ 2 \z) 
satisfy 


rf2 f (i,2) 

dz 2J+ 

p{z) 


pw|/ 1 1 ' 2, + 9 W/< 1 ' 2 » = o. 

1-gp l-g f 1 

2 2 — t Z — A 7 


?(*) 


1 Cq Ct // 

- -I—- -I- - —I— t- —, 

4 2 2 — £ 2 — A 1 


(29) 


where Co, C*, A and p are complicated functions of the 
entries of A(z). This EDO has, along with the singular 
points at 2 = 0,£, 00 , an extra singularity at 2 = A. This 
singularity can be checked to be an apparent one: the 
solutions of the indicial equation at 2 = A gives aj = 0, 2 
and there is no logarithm behavior due to the algebraic 


The tau-function has the direct interpretation of gener¬ 
ating function for correlations in field-theoretic applica¬ 
tions of the Painleve transcendents. Asymptotic expres¬ 
sions for the tau-function have been derived in 0 and 
the (irregular) conformal block interpretation was given 
in [23|, [24[ , and the relevant results are listed in the Ap¬ 
pendix. In order to describe it we define the composite 
monodromy parameter 

2COS7R7 = Tr(M* M 0 ) =Ti{M~ 1 ) 

= 2 cos 7T0OO + e 7 ”® 00 si s 2 , 

then the monodromy data can be written as 0 = 
{# 0 , 9 t , Oqo, ct, s^. We will assume generic (i.e., non¬ 
multiples of 7 r) values for the monodromy data 9 so these 
expressions can be locally inverted. 

In terms of (131 [I , the existence of monodromy matrices 
for the Rabi Fuchsian system m amounts to the exis¬ 
tence of a solution to the tau-function given the initial 
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conditions: 


— log r(t,6>) 


t=—4g 2 


E + g 2 A 2 

2 + V 


W ,ogT (t ’ S) 


t=-ig'■ 


1 A 2 

= -TT Tr A()A t = - 7 , 

t 2 16 g 4 ’ 


(33) 


which is guaranteed on general grounds. Note that in 
the case of interest 9 0 = 9 t = E + g 2 and 9^ = 0, these 
conditions can be inverted to yield the non-trivial mon- 
odromy data - the Stokes parameters in our application. 
It should be stressed that the equations above solve the 
Rabi model in an implicit but combinatorial sense: the 
formulae given in the Appendix give an asymptotic ex¬ 
pansion for the Painleve V tau-function near t = 0. The 
same special type of Painleve V system as above was 
studied in a series of papers 2^- 271, where the invariants 
of the isomonodromy flow were calculated, and, in the 
last of the series, a Toda chain structure outlined. 


QUANTIZATION 

The pair of equations (1331) provides an implicit solu¬ 
tion of the Rabi model in terms of the Painleve V tau- 
function. The quantization condition for the energy can 
be solved in a similar form by making use of the compos¬ 
ite monodromy parameter (l32l) . We start by considering 
the asymptotics of the solution of the Fuchsian system ([9]) 
with the behavior near z = oo given by (USD- Near the 
regular singular points z = 0, t the solution will behave 
as: 


ing that the matrix C$C t -1 which connects the natural 
solutions at z = 0 and z = t is diagonal: the analytic 
solution at 2 = 0 will also be analytic at z = t. In prin¬ 
ciple the connection could be “upper triangular”: the 
second solution at z = 0, which diverges as z B °, could 
be connected to a superposition of the divergent and the 
regular solutions at z = t, but one can easily see that 
this does not happen: consider the determinant of the 
fundamental matrix, det<f>, which satisfies the equation 

-j- det <&=(—-!—) det $. (35) 

dz \ z z — t) 

This equation yields det $ = z e °(z — t) St , and can be 
used to “normalize” the solutions, in the sense that now 
the fundamental matrix has unit determinant. One can 
convince oneself that this normalization does not change 
the connection matrices Gi, but now the two natural so¬ 
lutions at any particular singular points have a similar 
behavior: (z — Zi ) ±0i / 2 . This parity, a Z 2 parity, was 
fundamental in Braak’s work. In [22j, it was though as¬ 
sociated to a time-reversal symmetry. For our applica¬ 
tion, this symmetry guarantees that the vanishing of one 
of the off-diagonal elements of CqC f _1 will imply the van¬ 
ishing of the other off-diagonal term. Hence GoG t _1 will 
be diagonal. 

Now, a diagonal CqC f _1 implies, for the composite 
monodromy parameter < 7 , defined in (1321) . that 

cos7rcr = cos7r(0 o + 9 t ). (36) 

Therefore, the regularity of the solution is expressed as 
a quantization condition: 


(Goz^o o)(t + 0(z))C 0 , 0 

$(z) = { 

1 ( e t 0 \ 

[g*(z — f)( 0 o/(i -|- 0(z — t))Ct, z —> t, 

where Gi (i = 0,t) are the matrices diagonalizing Ai , 
and Ci are called the connection matrices. The mon¬ 
odromy matrices are diagonalized by the connection ma¬ 
trices, that is, 

M i = C~ 1 e niei C i . (34) 


(T n — 2n 0q 6t — 2(_£/ H- g ^ H- 77 .), n G Zu (37) 

Since a is given in terms of the Stokes parameters Si, 2 , 
this condition can be fed into the solution (1531) to yield 
the quantized values for the energy E n . The completion 
of this task requires the knowledge of the expansion of 
the tau-function given in the Appendix. 

PERSPECTIVES 


Thus, the Cj are defined up to right multiplication. With¬ 
out loss of generality we can take det Gj = 1. The matrix 
GoG^ 1 can be seen to “connect” the natural solutions of 
the system © at z = 0 and z = t. 

We can use the monodromy matrix to solve for the 
eigenvalue problem. The solution of © is required from 
physical grounds to be analytic on the whole plane - it 
will have an essential singularity at z = oo, but analytic- 
ity at z = 0, t ensures that the quantum state defined by 
the solution has finite expectation values for the relevant 
physical quantities (like the bosonic number operator). 
This condition is translated to our language by requir¬ 


The methods described here are useful not only to show 
the existence of the monodromy matrices, and hence 
Yang-Baxter integrability, for the Rabi model but it also 
provides a solution for the eigenvalue problem in terms of 
the transcendental equation (l33l) . Given that there is a 
combinatorial expansion of the Painleve V tau-function 
in terms of irregular conformal blocks, one can then im¬ 
plement a numerical/symbolic computation to complete 
the task of finding the eigenvalues, using the expansion 
given in the Appendix and the quantization condition 

(El). 

Another interesting direction would be to use the pro- 
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posed formalism to other similar systems. For instance, 
the extension to the model with broken parity introduced 
by Braak. This consists in adding a term of the form r )o x 
to the Rabi Hamiltonian, which seems to be a simple 
extension and amenable through the methods described 
here. 
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Formulae for the Painleve V r-function 


Here we lift the relevant formulae from [24]. In the 
following we consider the tau-function as defined in 53 : 




(38) 


The expansion for the tau-function is of the form: 
f(t,9)= ^2C({9i},a + n)s n t^ +n) B({9i}, a + n;t), 

n(=Z 

(39) 

where the irregular conformal block B is given as a power 
series over the set of Young tableaux ¥: 

B{{9 l } 1 a-t) = e "*** £ S Ai/i ({0 i }, t r)i |A|+ W (40) 

a,^gy 


p> _ tt (flop + a + i ~ + & + i ~ j) 2 - Aq) 

^ ( AA A A' +^-i-j + l + 2a) 

tt (^oo — cr + i - j)((9t - & + i ~ j) 2 - 9 q) 

hl(i > j)(X i +^ j -i-j + l + 2a) ' 1 j 

where A denotes a Young tableau, A* is the number of 
boxes in row i, A' is the number of boxes in column j 
and h\(i,j) = Aj + A'- — i — j + 1 is the hook length 
related to the box ( i,j ) G A. The structure constants G 
are rational products of Barnes functions: 

G{{9i\,cr) = G(1 + 0oo + e<r)G(l + 0* + 9q + ecr)x 

e=± 

G(1 + 9 t — 9 0 + ea)/G(l + 2ecr), (42) 

where G(z) is defined by the functional equation G(1 + 
z) = T(z)G(z). The parameters a and s in (l39l) are re¬ 


lated to the “constants of integration” of the Painleve V 
equation. The a is the same monodromy parameter as 
(1521) . whereas s has a rather lengthy expression in terms 
of monodromy data that can be read from 0- The 
Painleve V tau-function was also considered in great de¬ 
tail in [13]. The particular set of parameters considered 


here were also considered in 25-27]. 
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